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ANALYSIS OF THE FLIGHT MOTIONS OF A SMALL 
DEPLOYABLE GLIDER CONFIGURATION (U) 

COORD NO. AF -AM -4 19 

Paul L. Coe, Jr. 

Langley Research Center 

SUMMARY 

An investigation was conducted at the request of the U.S. Air Force Avionics 
Laboratory to determine the cause of a series of unsuccessful test flights of a small 
uncontrolled glider with folding wings. Initially, contractor tests conducted in a swim- 
ming pool indicated that the glider exhibited satisfactory glide capability, trim charac- 
teristics, and dynamic stability characteristics. The pool tests were followed by a num- 
ber of unsuccessful flight tests in air during which the model generally exhibited tight 
spirals or continuous 360° rolling motions which resulted in an essentially ballistic 
trajectory. 

The present study consisted of wind-tunnel tests of an actual glider and a theoretical 
analysis of the performance, stability, and trimmability of the configuration. The analysis 
also defined the factors which probably caused the vastly different motions exhibited by the 
configuration in air and in water. 

The results of the study showed that the proximity of the trim condition to the onset 
of wing stall, the inherent lateral-directional instability of the configuration, and the exist- 
ence of large out-of-trim rolling moments created by wing asymmetries were all possible 
causes for the series of unsuccessful test flights. 

INTRODUCTION 

An investigation was conducted at the request of the U.S. Air Force Avionics 
Laboratory to determine the cause of a series of unsuccessful test flights of a small 
uncontrolled glider with folding wings. In operational use, it was intended that a large 
number of the gliders, containing electronic countermeasure (ECM) equipment, would be 
simultaneously deployed from aircraft during tactical missions. Initially, contractor 
tests were conducted by launching the vehicle under water in a swimming pool in order 
to determine the performance, and stability of the glider motions. The results of these 
pool tests indicated that the glider exhibited satisfactory glide capability, trim character- 
istics, and dynamic stability characteristics. The pool tests were followed by a number 
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of unsuccessful flight tests in air involving drop tests from a general-aviation airplane 
at an altitude of 609.6 m (2000 ft). After the launch, the model generally exhibited tight 
spirals or continuous 360° rolling motions which resulted in an essentially ballistic tra- 
jectory. In some instances, however, a few of the gliders appeared to recover from the 
launch and obtain a momentary wings-level glide that was followed immediately by a sud- 
den wing drop and 360° rolls along a ballistic trajectory. In addition, several flight tests 
were made with a model with a smaller vertical tail. In these tests the model generally 
exhibited the same rolling motions; however, a few of the models were seen to exhibit a 
severe Dutch roll oscillation. As a result of the rapid rolling motions, the various gliders 
tested exhibited essentially no glide capability, and the configuration was obviously unsuit- 
able for the intended operational mission. 

The present study, conducted at the NASA Langley Research Center, consisted of 
wind-tunnel tests of an actual glider and a theoretical analysis of the performance, sta- 
bility, and trimmability of the configuration. The analysis also defined the factors which 
probably caused the vastly different motions exhibited by the gliders in air and in water. 

SYMBOLS 

All aerodynamic data with the exception of lift and drag are presented with respect 
to the body system of axes. Moment data are presented with respect to a center-of-gravity 
position of 25 percent of the wing mean aerodynamic chord. Dimensional values are pre- 
sented in the International System of Units (SI) with equivalent values given parenthetically 
in the U.S. Customary Units. 

b wing span, m (ft) 

Cj) drag coefficient, Fj^/q^S 

Cl lift coefficient, Fl/Q^S 

C^ rolling-moment coefficient, Mx/^^Sb 

o rolling- moment coefficient at /3 = 0° 

C m pitching-moment coefficient, My/q^Sc 

C n yawing-moment coefficient, Mg/q^Sb 

Cx longitudinal-force coefficient, Fx/q^S 
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side-force coefficient, Fy/qJS 

normal-force coefficient, Fg/q^S 

wing mean aerodynamic chord, m (ft) 

drag force, N (lbf) 

lift force, N (lbf) 

longitudinal force, N (lbf) 

side force, N (lbf) 

normal force, N (lbf) 

altitude, m (ft) 

moment of inertia about longitudinal body axis, kg 
moment of inertia about lateral body axis, kg-m2 
moment of inertia about normal body axis, kg-m2 
lift- drag ratio 

rolling moment, m-N (ft-lbf) 

pitching moment, m-N (ft-lbf) 

yawing moment, m-N (ft-lbf) 

glider mass, kg (slugs) 

rolling velocity, rad/sec or deg/sec 

pitching velocity, rad/sec or deg/sec 

free- stream dynamic pressure, N/m2 (lbf/ft 2 ) 



m2 (slug-ft2) 

(slug-ft2) 

(slug-ft2) 


r 


yawing velocity, rad/sec or deg/sec 


S wing area, m 2 (ft 2 ) 

t time, sec 

tj/2 time required for amplitude of oscillation to decrease by a factor of 2, sec 

V velocity, m/sec (ft/sec) 

o? angle of attack, deg or rad 

angle of sideslip, deg or rad 
y flight-path angle, deg or rad 

0 inclination of X body axis with respect to horizon 

M relative density parameter, m/pSb 

p mass density of air or water, kg/m 3 (slugs/ft 3 ) 

o angle of bank, deg or rad 

ip angle of yaw, deg or rad 

Coefficients: 


'Xq 


Cl fi dft 


aCx 

c Zq = 

aCz 

Cm a - 

aC m 

qc 

a — 
2V 


qc 

a — 
2V 

q 

qc 
9 — 
2V 

8C Z 

c no = 

ac n 
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CY fl = 

bc y 

r\ A 


dC 1 

p EeE 

2V 


C *r = 


Cl 
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DESCRIPTION OF GLIDER 


A three -view sketch showing the general layout of the glider configuration is pre- 
sented in figure 1 and the mass and geometric characteristics of the vehicle are given in 
table I. The configuration utilized deployable (spring-loaded) aluminum sheet metal wings, 
a hardwood fuselage, and sheet metal tail surfaces. As indicated in table I, the nominal 
glider mass of 0.68 kg (0.0466 slug) resulted in a high value of the relative density p 
(approximately 93.3 for the vehicle in air at 609.6 m (2000 ft)). 

One of the most important geometric features of the glider was the wing construction 
and shape. The wing panels were stamped aluminum slabs shaped to a 10-percent mean 
camber line. The techniques involved in the construction of the wings and mating of the 
wing panels to the fuselage were not under precise quality control, and it will subsequently 
be shown that aerodynamic characteristics arising from model asymmetries were an 
important factor in the unsuccessful flight tests. 

METHOD OF ANALYSIS 

Static wind-tunnel force tests were conducted to determine the longitudinal and lat- 
eral directional aerodynamic characteristics of the glider for the two vertical-tail config- 
urations shown in figure 1. Dynamic stability derivatives of the glider were estimated by 
use of the methods of reference 1. The static and dynamic aerodynamic stability deriva- 
tives were used as input data for theoretical calculations of the lateral-directional dynamic 
stability characteristics using linear three-degrees-of-freedom equations of motion. In 
addition, time histories of flight motions were generated by use of a nonlinear six-degrees- 
uf-fri i dun. digital computer program. 


Wind Tunnel 

The static aerodynamic data presented herein were obtained at the Langley Research 
Center with a representative glider in a low-speed wind tunnel with a 3.66-m (12-ft) octag- 
onal test section. The model was strut mounted to a small six-component strain-gage bal- 
ance. Tests were conducted with the glider erect and inverted in order to determine longi- 
tudinal flow angularity, and a flow survey was made to determine the directional flow 
angularity. The test apparatus was alined with the flow so that the glider was at ( 3 = 0 ° 
for 0 = 0°. 


Tests 

Static force tests were conducted in the low-speed wind tunnel at a Reynolds number 
of approximately 0.04 x 10 6 based on the mean aerodynamic chord of the wing. During the 
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tests, measurements were made of the six force and moment components for a center-of- 
gravity location corresponding to 0.2 5c over an angle-of-attack range of -10 to 20 for a 
range of angles of sideslip of ±10°. 


Calculations 

The dynamic stability characteristics of the lateral-directional modes of motion 
were calculated for the glider with each of the vertical tails shown in figure 1 by using 
the linearized three -degrees-of -freedom equations of motion presented in reference 2. 

The calculated characteristics included the time to half-amplitude t\/2 the Dutch 
roll mode, spiral mode, and roll-subsidence modes. 

Motions of the glider in air and in water were calculated by using a nonlinear, six- 
degrees-of-f reedom computer program. Use of the nonlinear program was necessitated 
by the large, rapid angular motions exhibited by the glider in air. 

RESULTS OF FORCE TESTS 
Longitudinal Characteristics 

The static longitudinal aerodynamic characteristics of the configuration are pre- 
sented in figure 2 for the complete configuration and the configuration with the horizontal 
tails removed. The data show that the basic configuration exhibited good longitudinal 
stability (negative, linear variation of Cm with a) and that the horizontal tail provided 
a contribution to longitudinal stability which was essentially constant over the angle-of- 
attack range of the tests. The data also indicate that the glider would trim at an angle of 
attack of about 2° with a value of L/D of about 5.2. During the longitudinal-force tests, 
a limited number of tuft flow-visualization tests were conducted to determine the flow con- 
ditions over the wing panels. The results of these flow -visualization tests indicated that 
the onset of wing stall occurred near a = 2°; this result is in agreement with the static 
longitudinal data presented in figure 2. The flow -visualization tests also indicated that 
the wings were completely stalled at a = 5°. Thus, the proximity of the trim angle of 
attack to the stall may have been the factor which resulted in the sudden wing drop follow- 
ing an initial pullout from launch. 

Lateral-Directional Characteristics 

Some of the more important lateral-directional aerodynamic characteristics of the 
glider are presented in figure 3 for the configuration with the small and large vertical 
tails and for the configuration with the vertical tail removed. The data are presented in 
terms of the static directional stability derivative C n ^ and the effective dihedral deriv- 
ative C ip, obtained over a range of sideslip angle of ±5°. Figure 3 also presents the 
variation of rolling- moment coefficient at zero sideslip C^q with angle of attack. 
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The variations of C n ^ and C ^ with a show that the configuration exhibited 
positive directional stability ^positive values of C n ^) and positive effective dihedral 
(negative values of over the range of angle of attack with either vertical tail. 


The factor most nearly relating to observed flight behavior difficulties was the 
indication that significant out-of-trim rolling moments existed for the configuration 
within the intended operational range of angle of attack. As shown in figure 3, the mag- 
nitude of the rolling-moment coefficient at /3 = 0° was large over the region of interest. 
For purposes of comparison, the magnitude of the rolling-moment coefficient measured 
at « = 0° was about equal to the magnitude obtained for full aileron deflection on con- 
ventional airplanes. The large out-of-trim values of C^q were probably caused by 
difficulties associated with holding close tolerances by use of simple manufacturing 
techniques. 


RESULTS OF CALCULATIONS 
Dynamic Stability Characteristics 

The results of calculations to determine the lateral-directional dynamic stability 
characteristics of the configuration for the two vertical-tail sizes used by the contractor 
are presented in figure 4. The experimentally determined values of the static derivatives, 
and the estimated values of the dynamic derivatives (see ref. 1) used in the calculations of 
the dynamic stability characteristics herein, are presented in table II. 

The results of the calculations for trimmed flight at various angles of attack for 
an altitude of 609.6 m (2000 ft) are presented in figure 4 in terms of the damping factor 
l/ti /2 of the roll, spiral, and Dutch roll modes. Positive values of l/tj /2 represent 
damped (dynamically stable) modes whereas negative values represent undamped (dynam- 
ically unstable) modes. The results of these calculations show that with the large verti- 
cal tail (fig. 4(a)), the configuration would be expected to exhibit an unstable spiral mode 
throughout the angle -of-attack range considered, and that the stability of all the modes 
decreased with increasing angle of attack. With the small vertical tail (fig. 4(b)), the 
results indicate the existence of an unstable Dutch roll mode for angles of attack from 
0° to 6°. An additional factor relating to the unsuccessful free-f light tests therefore 
appears to be the inherent dynamic instability of the configuration. 

Presented in figure 5 are boundaries for spiral stability and oscillatory (Dutch roll) 
stability in terms of variations of the stability derivatives C n ? and for the config- 

uration at oi — 2°. Values of Cn.^ and C ^ are also indicated for the glider with the 
small and large vertical tails. The locations of these values relative to the spiral and 
oscillatory boundaries further illustrate the inherent dynamic stability problem, and 
the data also indicate that a stable region can be obtained with suitable values of C n ^ 
and Cip. 
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The data of figure 5 also indicate that the dynamic stability characteristics exhibited 
by the configuration may be vastly different in water. Although the spiral stability bound- 
ary is independent of the relative density of the glider, the dashed line in figure 5 shows 
that the oscillatory stability boundary becomes markedly less restrictive, with the result 
that there is a much larger region of stability. It would appear, therefore, that the glider 
configuration would exhibit unrealistically good oscillatory stability in the swimming pool 
tests. 


Calculated Flight Motions 

Presented in figures 6,7, and 8 are calculated time histories of flight motions for 
the glider in air and in water. The initial conditions for the calculations are as follows: 


Condition 

Initial condition in — 

Air 

Water 

h, m (ft) 

609.6 (2000) 

1.83 (6) 

V, m/sec (ft/sec) 

30.48 (100) 

0.91 (3) 

a, deg 

3 

3 

ft, deg 

1 

1 

p, kg/m3 (slugs/ft3) .... 

1.154 (0.00224) 

999.66 (1.94) 

0, deg 

0 

0 

0, deg 

-3 

-3 

0, deg 

0 

i 

0 


Shown in figure 6(a) are the calculated motions for the configuration with the large 
vertical tail. The Dutch roll oscillation is seen to be stable (as expected based on the 
data of fig. 4) as evidenced by the variations in the heading i//. After the initial distur- 
bance, the motions damp out and the glider trims at a - 2°. As time progresses, the 
spiral instability previously discussed for this configuration begins to appear in the form 
of continuously increasing values of bank angle 0 and heading 0. The total flight time 
from launch at h = 609.6 m (2000 ft) to impact at h = 0 is approximately 142 seconds. 

The time history for the glider with the small vertical tail is presented in figure 6(b). 
The motion is characterized initially by the extremely unstable Dutch roll oscillation, pre- 
viously discussed for this configuration, which is seen to build up to a limit cycle oscilla- 
tion for which variations in bank angles of ±80° occurred. This condition resulted in a 
significant decrease in total flight time to approximately 78 seconds (as compared with 
142 seconds for the configuration with the larger vertical tail). 
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Although the foregoing calculated motions may appear to be acceptable for the 
glider, it must be remembered that the effects of asymmetric wing stall have not been 
simulated. Additionally, these results have been obtained by use of relatively small 
disturbances from the trim condition, whereas larger disturbances would be expected 
to yield significantly larger amplitude motions which may result in completely unaccept- 
able behavior. 

Figure 7 presents the results of time history calculations to determine the effects 
of out-of-trim rolling moments. The calculations were made for the configuration with 
both the large and small vertical tails and included the values of out-of-trim rolling 
moment shown in figure 3. The results of the calculations show that for both tail sizes 
considered, the configuration immediately developed a large roll rate which resulted in 
a continuous 360° rolling motion. As a result of the rolling motion, the glider followed 
an essentially ballistic trajectory and reached h = 0 at t = 14.5 sec. This result is in 
good agreement with flight-test results and appears to indicate that the primary factor 
causing the rapid rolling motion was the out-of-trim rolling moment. 

Since the models with the small vertical tail were observed to exhibit satisfactory 
behavior when tested in water, additional time histories were computed to simulate this 
condition. The time histories in figure 8(a) show the motion of the glider in water for 
the case of no rolling-moment asymmetry. The configuration exhibits excellent flight 
characteristics with a heavily damped Dutch roll oscillation, as expected, based on the 
data of figure 5 which show the marked shift in the oscillatory stability boundary when 
the density increases from that of air to that of water. Figure 8(b) shows the L ime his- 
tories for the configuration in water with the out-of-trim rolling moments included. The 
configuration developed a relatively small roll rate and exhibited steady turning flight 
with a bank angle of approximately 1°. 

The differences between the mild turn exhibited by the configuration in a ter and 
the ballistic flight path determined in air can be explained by consideration of the rela- 
tionship between the rolling-moment coefficient required for steady turning flight C^q 
and the relative density parameter p . An approximate expression showing this relation- 
ship is developed in the appendix as 



Using the measured and estimated values of the stability derivatives at the trim con- 
dition thus defines this relationship which is plotted in figure 9 with a parametr c depend- 
ence on the bank angle 0. From figure 9 it is seen that for the level of asymr 'try con- 
sidered herein, the vehicle would be expected to turn with a bank angle of apprc mately 1° 
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in water. This result correlates with the time history calculations shown in figure 8(b). 
Additionally, from figure 9 the vastly larger relative density of the vehicle in air places 
the configuration in a region for which no steady- state equilibrium turn is possible. Thus, 
in air, the out-of-trim rolling moment would be expected to cause the vehicle to develop a 
large body axis roll rate, which would then result in an essentially ballistic trajectory. 
This result correlates with the time history calculations shown in figure 7, and with the 
experimental flight-test results. 

SUMMARY OF RESULTS 

The results of an investigation to determine the stability characteristics of a deploy- 
able small glider, with folding wings, have indicated the following possible causes for a 
series of unsuccessful test flights: 

1. The proximity of the trim condition to the onset of wing stall 

2. The inherent lateral-directional instability of the configuration 

3. The existence of large out-of-trim rolling moments created by wing asymmetries. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 28, 1975. 
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APPENDIX 


ROLLING MOMENT REQUIRED FOR STEADY TURNING FLIGHT 


An approximate expression for the rolling- moment coefficient required for steady 
turning flight C^q as a function of the relative density parameter p is developed 
through consideration of the forces and moments acting on the vehicle shown in figure 10 
where R is the radius of turn, y is the flight-path angle, and J2 is the turn rate. 
Centripetal equilibrium requires 

Fl cos (-y) + Fj) sin (-y)J sin 0 + Fy cos 0 = - (Al) 

Nondimensionalizing equation (Al) and neglecting the side force yields 


C L cos (-y) + Cd sin (-y) 


sin 0 


2m 

RpS 


(A2) 


An approximate form of tlie yawing- and rolling-moment equations may be written as 
follows: 


C„^ + C„ p j| + C„ r £|*0 

- C ip § + C lr § + C i>0 - 0 


If small values of 6 are assumed, the following approximate expressions for the non- 
dimensional rolling and yawing velocities result: 

EL _ . £b (A5) 

2V 2V 


rb 

2V 


fib , 
2V C0S * 


The assumption of small values of y yields 


(AC; 


_ _b_ (AT) 

2V 2R 

Substituting equation (A7) into equation (A2) and equations (A5) and (A6) into equations (A3) 
and (A4) yields the following set of approximate expressions: 

[Cl cos (-y) + Cd sin (-y)j sin 0 ~ 4 p ( A8 ) 
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APPENDIX - Concluded 

CnpP + (c nr cos 0 - C np e)p£ * 0 (A9) 


Cipft + (c lr cos 0 - C lp 0 )|| + C lf0 « 0 (A10) 

Solving equations (A9) and (A10) for fJb/2V and substituting this resu t into equation (A8) 
yields upon rearranging 


] V -L 


C 


*,0 " 


Cl{fn r ~ c n j 3 C Z r ) cos ~ c lp c n { 


sin 0 


4p C n , 


(All) 


Equation (All) is an approximate expression for the rolling-moment coefficient required 
for steady turning flight C^q as a function of the relative density parameter p and the 
bank angle 0. 
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TABLE I.- MASS AND GEOMETRIC CHARACTERISTICS OF THE GLIDER 


Mass (nominal), kg (slugs) 

Moments of inertia, kg~m2 (slug-ft2) 

l X 

IY 

iz 

Overall length, cm (in.) 

Wing: 

Span, cm (in.) 

Area, cm^ (in^) 

Mean aerodynamic chord, cm (in.) 

Aspect ratio 

Taper ratio 

Dihedral, deg 

Incidence, deg 

Airfoil section 

Horizontal tail: 

Span, cm (in.) 

Area, cm2 (in2) 

Mean aerodynamic chord, cm (in.) 

Aspect ratio 

Taper ratio 

Dihedral 

Incidence, deg 

Airfoil section 

Large vertical tail: 

Span, cm (in.) 

Area, cra^ (in^) 

Root chord, cm (in.) 

Tip chord, cm (in.) 

Airfoil section 


0.68 (0.0466) 

0.0020 (0.0015) 

0.0027 (0.0020) 

0.0045 (0.0033) 

26.16 (10.3) 

40.64 (16.0) 

154.84 (24.0) 

3.81 (1.5) 

10.67 

1.0 

3.0 

3.5 

10-percent-camber line metal slab 


15.24 (6.0) 

58.06 (9.0) 
3.81 (1.5) 
. . . . 4.0 

. . . . 1.0 

. . . . 0 

. . . . -1.5 

. Metal slab 

6.35 (2.5) 
26.19 (4.06 '• 
4.70 (1.85) 
3.56 (1.4) 
Metal slab 


Small vertical tail: 

Span, cm (in.) 3.81 (1.5) 

Area, cm2 (i n 2) 16.19 (2.51) 

Root chord, cm (in.) 4.70 (1.85) 

Tip chord, cm (in.) 3.81 (1.5) 

Airfoil section Metal slab 
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TABLE II.- STABILITY DERIVATIVES 


(a) Lateral 


Derivative 

Value for — 

Large vertical tail at o 

of - 

Small vertical tail at a 

of - 

0° 

4° 

8° 

12° 

0° 

4° 

8° 

12° 

C Y0 

-0.470 

-0.470 

-0.540 

-0.610 

-0.068 

-0.068 

-0.202 

-0.270 

Cn 

.200 

.200 

.230 

.260 

.029 

.029 

.086 

.115 

% 

-.172 

-.200 

-.172 

-.200 

-.200 

-.200 

-.212 

-.218 

C Yp 

0 

0 

0 


0 

0 

0 

0 

Cn p 

-.05 

0 


.06 

-.05 

0 

.06 

.06 

^Zp 

-.70 

-.37 

-.26 

-.06 

-.70 

-.37 

-.26 

-.06 

CY r 

.40 

.40 

.46 

.52 

.06 

.06 

.17 

.23 

Cn r 

-.28 

-.29 

-.33 

-.36 

-.14 

-.14 

-.20 

-.23 

c £r 

.25 

.40 

.43 

.47 

.27 

.40 

.48 

.49 


C 




= 0 


(b) Longitudinal 
c Z q = ° 


C mq — 45.0 
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Figure 5.- Lateral stability boundaries. 




? r JIJIV 



r 5 r 

Z l/W-— 
oL 


s 



-bL 



5^ 



ior 

? 5 — 

oiA/Vw^ 



hV 

- I 2S00 I [ 

oL L. _J_ __1 _1 ..._L J .i i L L I I I 1 1 

0 2 4 B 8 10 12 14 16 18 20 22 24 26 28 30 

T i -e , set 

(a) Configuration with large vertical tail. 

Figure 6.- Calculated time history of flight motions in air. (Out-of-trim rolling 

moment not included.) 
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(b) Configuration with small vertical tail. 
Figure 6.- Concluded. 
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(a) Configuration with large vertical tail. 

Figure 7.- Calculated time history of flight motions in air. (Out-of-trim rolling 

moment included.) 
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(a) Out-of-trim rolling moment not included. 

Figure 8.- Calculated time history of motions in water. (Configuration with 

small vertical tail.) 
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(b) Out-of-trim rolling moment included. 


Figure 8.- Concluded. 
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Figure 9.- Approximate variation of rolling -moment coefficient required for steady -state 
turn. Configuration with small vertical tail; a = 2°. 
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